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1.0 


INTRODUCTION 


1 . 1 Background 


This  memorandum  documents  recent  improvements  to  t+rg  model 


developed  by  Bolt  Beranek  and  Newman  (BBN)vfor  the  prediction  of 
volume-time  histories  for  sheet  cavitation  on  the  blades  of  mer- 
chant ship  propellers,  and  the  low  frequency  acoustic  radiation 
resulting  from  these  sheet  cavitation  volume  fluctuations.  This 
model  is  part  of  a larger  effort  directed  at  predicting  the  stat- 
istics of  ocean  ambient  noise  due  to  merchant  shipping. 


A description  of  the  mechanics  of  sheet  cavitation  on  pro- 
peller blades  and  the  acoustic  radiation  resulting  from  fluctua- 
tions in  sheet  cavitation  volume  is  given  in  an  earlier  report, 
BBN  Tech  Memo  No.  319,  and  will  not  be  repeated  herev 

It  is  assumed  that  the  reader  is  familiar  with  the  contents  of 
Tech  Memo  No.  319,  which  describes  these  processes  and  the  ap- 
proach used  to  model  them. 


1.2  Objectives  and  Scope 


The  previous  version  of  our  model  considerably  over-predicted 
the  amount  of  cavitation  present  on  propeller  blades  as  they  ro- 
tate through  the  uneven  wakefield  behind  a single  screw  ship,  and 
thus  over-predicted  the  acoustic  radiation  resulting  from  these 
cavity  fluctuations.  As  discussed  in  Reference  2,  this  over-pre- 
diction  was  attributed  to  three  major  factors: 

o For  determining  the  loading  (lift)  on  the  propeller  blade, 
the  previous  model  assumed  that  the  blade  could  be  replaced 
by  a single  "lifting  line"  of  zero  chord,  and  that  the 
inflow  at  this  lifting  line  was  representative  of  the  flow 
over  the  whole  blade.  Physical  propellers  have  finite 
chord  lengths  and  tend  to  "average"  the  inflow  variations 
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over  the  chordwise  extent  of  the  blade,  thus  reducing 
the  response  of  the  blade  to  inflow  variations. 

o The  previous  model  assumed  that  the  propeller  blade  sec- 
tion could  be  replaced  by  a zero  thickness  flat  plate 
when  calculating  the  occurrence  and  extent  of  shett  cavi- 
tation. This  will  tend  to  over-predict  the  occurrence  of 
cavitation,  since  real  propellers  are  comprised  of  airfoil 
type  sections  with  rounded  leading  edges  and  can  tolerate 
moderate  angle  of  attack  fluctuations  without  provoking 
the  occurrence  of  sheet  cavitation. 

o The  previous  model  did  not  include  the  effect  of  cavita- 
tion on  the  lift  of  the  propeller  blade  sections  (cavitat- 
ing  lift-slope ) . 

The  present  report  documents  four  major  items: 

o Improvements  to  the  model  to  overcome  the  shortcomings 
listed  above. 

o A sensitivity  analysis  to  determine  the  influence  of  the 
input  variables  on  the  predicted  cavitation  noise. 

o Comparison  of  measured  cavitation  line  spectra  from  mer- 
chant ships  with  levels  estimated  using  the  current  model. 

o Discussion  of  shortcomings  in  the  modeling  process  which 
require  further  investigation. 
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2.0  TECHNICAL  APPROACH 

The  approach  adopted  for  the  calculation  of  the  blade  loading 
at  each  angular  position  in  the  wake  is  based  on  quasi-steady  lift- 
ing-line theory,  as  described  in  Ref.  1.  The  details  of  the  numer- 
ical lifting-line  calculations  have  been  changed  somewhat  from  the 
description  given  in  Ref.  1 due  to  numerical  problems  in  the  pre- 
vious computer  program.  The  new  program  has  proven  to  be  very  re- 
liable and  shows  rapid  convergence  to  the  solution  of  the  inverse 
lifting  line  problem  with  no  numerical  instabilities  evident. 

2.1  Effects  of  Finite  Blade  Width 

In  order  to  account  for  the  effect  of  finite  blade  width 
when  calculating  the  loading  on  the  blade,  the  inflow  velocities 
are  averaged  over  the  chordwise  extent  of  the  blade  in  order  to 
determine  the  axial  wake  velocity  V^,  the  tangential  wake  velocity 
VT,  and  the  geometric  Inflow  angle  g (see  Fig.  2.1).  Initially 
the  averaging  was  done  uniformly  over  the  blade  chord,  i.e.  the 
wake  velocities  near  the  leading  and  trailing  edges  of  the  blade 
were  weighted  equally  in  the  averaging  process.  This  did  not  seem 
to  provide  the  correct  response  of  the  propeller  blade  to  the  wake- 
field.  Since  Brown  (3)  had  shown  the  leading  edge  of  the  propeller 
blade  to  govern  the  generation  of  unsteady  propeller  forces,  the 
averaging  calculations  were  modified  to  weight  the  inflow  veloci- 
ties near  the  leading  edge  of  the  blade  more  heavily  than  those 
near  the  trailing  edge.  This  appears  to  provide  a better  estimate 
of  the  propeller  blade's  response  to  the  wake  non-uniformities. 
Sasajima  (4)  has  also  applied  this  leading-edge  weighting  method 
to  the  quasi-steady  estimation  of  unsteady  propeller  forces  and 
found  generally  good  agreement  with  experimental  results  and  cal- 
culations done  using  very  complete  unsteady  lifting  surface  theory 
(5).  This  wake  averaging  method  with  the  averaging  weighted  toward 
the  leading  edge  of  the  blade,  has  been  adopted  for  the  present  work. 
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2.2  Propeller  Blade  Sections 

Actual  propeller  blade  sections  are  able  to  withstand  a finite 
change  In  angle  of  attack  (Aa)  without  suffering  leading  edge  sheet 
cavitation.  The  angle  of  attack,  a that  a given  section  can  with- 
stand is  dependent  on  the  cavitation  number,  o 

P -P„ 

•'ipvi-  <2-1> 


p = 

O 

PV  = 
P = 

V = 


ambient  pressure 
vapor  pressure  of  fluid 
fluid  density 
flow  velocity 


Graphically,  this  is  represented  by  a cavitation  "bucket"  diagram, 
shown  in  Fig.  2.2.  The  "bucket"  is  the  curve  drawn  on  a a vs  a 
diagram  showing  the  boundary  between  the  region  for  which  sheet 
cavitation  will  occur  and  the  region  of  no  cavitation.  Note  that 
the  bucket  width  for  a given  section  is  dependent  on  the  cavitation 
number  o. 

The  shape  of  the  cavitation  bucket  for  a propeller  blade  sec- 
tion is  dependent  on  the  thickness-to-chord  (t/c)  ratio,  the  type 
of  thickness  distribution,  and  the  magnitude  and  shape  of  the  cam- 
ber distribution  (curvature  of  the  centerline  of  the  blade  section). 
For  a given  arbitrary  blade  section  the  bucket  curve  may  be  calcu- 
lated, but  the  numerical  calculations  are  quite  lengthy,  and  out- 
side the  scope  of  the  present  work.  The  approach  utilized  here 
is  to  consult  tabulated  cavitation  buckets  for  a given  type  of 
thickness  and  camber  distribution  and  assume  that  these  sections 
are  representative  of  the  type  of  blade  sections  found  in  the 
world's  merchant  ship  propellers.  The  type  of  section  used  in  the 
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present  work  is  the  NACA  66  thickness  distribution  in  conjunction 
with  an  a=0.8  type  of  camber  distribution  (Ref.  6).  This  type  of 
section  is  widely  used  in  marine  propellers  designed  and  built  in 
the  past  few  years.  Very  complete  cavitation  bucket  diagrams  are 
available  for  this  section  (Ref.  7),  and  the  buckets  were  already 
available  in  numerical  form  (8)  for  use  in  our  computer  program, 
which  saved  considerable  programming  effort  in  the  present  work. 

The  NACA  66  section,  whose  properties  have  been  incorporated 
into  the  present  computer  program,  is  more  tolerant  of  angle  of 
attack  fluctuations  than  many  types  of  blade  sections  used  in  the 
world's  merchant  ship  propellers.  However,  since  the  thickness- 
to-chord  (t/c)  ratio  is  the  major  parameter  which  affects  cavita- 
tion bucket  width  (Aa)  at  a given  cavitation  number  (a),  the  cavi- 
tation buckets  of  other  types  of  blade  sections  may  be  simulated 
by  entering  the  NACA  66  bucket  curves  at  a t/c  ratio  which  is 
suitabily  modified  so  that  the  resulting  bucket  width  is  that  of 
the  section  under  consideration.  This  modified  t/c  ratio  is 
readily  calculated. 

In  the  computer  program,  the  angle  of  attack  a for  each  blade 
section  of  the  propeller  is  determined  by  using  quasi-steady  lift- 
ing line  theory  to  solve  for  the  loading  on  the  blade,  as  discussed 
in  Sections  2.0  and  2.1.  The  cavitation  number  a for  each  section 
is  calculated  from  Equation  2.1.  Then  for  each  blade  section  along 
the  blade,  the  NACA  66  bucket  curves  are  entered  with  the  proper 
values  of  a,  a,  t/c,  and  camber  ratio  to  determine  whether  or  not 
sheet  cavitation  is  present  at  that  section. 

2.3  Calculation  of  Cavity  Extent  and  Cross  Sectional  Area 

Once  it  has  been  determined  that  sheet  cavitation  is  occurring 
at  a particular  blade  section,  the  chordwise  extent  and  cross-sec- 
tional area  of  the  cavity  at  that  radius  are  calculated  using  lin- 
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earized  free-streamline  cavitating  foil  theory  (9-13).  This 
theory  assumes  that  the  pressure  field  around  the  blade  section 
which  influences  the  cavitation  extent  and  area  is  composed  of 
two  components:  the  pressure  distribution  associated  with  a flat 
plate  at  angle  of  incidence  a,  and  the  pressure  distribution  assoc- 
iated with  a foil  having  a parabolic  arc  camberline  of  camber 
ratio  X,  as  shown  in  Figure  2.3.  The  effective  camber  of  a blade 
section  in  non-uniform  flow  can  be  considered  as  being  made  up 
of  two  parts: 

o The  geometric  camber  of  the  actual  blade  section. 

o The  induced  camber  caused  by  the  curvature  of  the  flow 
over  the  blade  chord. 

Van  Oossanen  has  shown  (14)  that  induced  camber  is  quite  import- 
ant in  determining  the  cavitation  characteristics  of  a propeller. 

We  have  adopted  van  Oossanen' s method  of  calculating  induced 
camber  due  to  non-uniform  inflow  (15),  and  use  the  sum  of  the  blade 
geometric  camber  and  the  induced  camber  in  entering  the  cavitating 
foil  calculations.  The  inclusion  of  this  induced  camber  effect 
significantly  improved  the  correlation  between  cavitation  patterns 
predicted  using  our  computer  program  and  the  cavitation  patterns 
observed  on  model  propellers  in  cavitation  tunnels. 

2.4  Effect  of  Cavitation  on  Lift  (Cavitating  Lift  Slope) 

Cavitation  is  known  to  alter  the  lift  and  drag  characteristics 
of  foil  sections  relative  to  the  non-cavitating  condition.  Since 
the  relationship  between  angle  of  attack  and  lift  coefficient  is 
used  to  solve  for  the  loading  on  the  blade,  the  effects  of  cavita- 
tion on  lift  slope  should  (in  principle)  be  included  in  the  quasi- 
steady solution  for  the  blade  loading.  We  originally  anticipated 
that  the  inclusion  of  this  effect  would  significantly  alter  the 
predicted  cavitation  volume  history  on  the  propeller.  A re-examin- 
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ation  of  this  problem  during  the  present  work  now  indicates  that 
the  effect  of  the  cavitating  lift  slope  on  the  final  predicted 
volume  history  would  be  quite  small,  and  the  error  involved  with 
neglecting  this  effect  is  Insignificant  compared  to  other  short- 
comings inherent  in  the  present  method  of  estimating  cavitation 
volume  histories.  Moreover,  the  inclusion  of  the  cavitating  lift 
slope  leads  to  numerical  instabilities  in  the  lifting-line  equa- 
tions which  may  be  very  difficult  to  overcome.  This  author  knows 
of  no  investigators  who  have  successfully  included  the  effects  of 
cavitation  in  solving  for  the  loading  on  the  propeller  blades. 

For  these  reasons  the  effect  of  the  cavitating  lift  slope  has  not 
been  included  in  the  present  work,  and  it  is  recommended  that  this 
phenomenon  be  ignored  until  other,  more  important  effects  have 
been  accounted  for  in  the  estimation  of  sheet  cavitation  volume- 
time histories. 

2.5  Summary  of  Analysis  Procedure 

The  procedure  used  to  calculate  the  cavitation  volume  on  a 
propeller  blade  at  a given  position  in  the  wakefield  is  summarized 
below . 

A.  For  each  of  9 blade  sections  along  the  blade,  examine  the 
wake  velocities  over  the  chordwise  extent  of  the  blade  and 
determine  the  weighted  average  inflow  angle  3,  axial  velo- 
city V^,  and  tangential  velocity  V,p.  The  induced  camber 
due  to  the  curvature  of  the  flow  is  also  calculated  at 
this  time. 

B.  Using  the  inflow  quantities  calculated  in  Step  A and  the 
known  blade  geometry,  calculate  the  blade  loading  and  cor- 
responding blade  section  angles  of  attack  using  quasi- 
steady lifting  line  theory. 
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C.  For  each  blade  section,  consult  the  appropriate  cavita- 
tion bucket  curve  to  determine  whether  or  not  sheet 
cavitation  is  occurring  at  that  section. 

D.  For  each  blade  section  at  which  sheet  cavitation  is  pre- 
dicted to  occurr,  calculate  the  chordwise  extent  and  the 
cross  sectional  area  of  the  sheet  cavity  using  free- 
streamline  cavitating  foil  theory. 

E.  Integrate  the  cavity  cross  sectional  areas  in  the  span- 
wise  direction  to  arrive  at  the  cavity  volume  on  the 
blade  for  a given  position  of  the  blade  in  the  wake. 

This  procedure  is  repeated  at  10°  increments  of  blade  rotation, 
so  that  the  volume  history  for  one  blade  may  be  determined  over 
360°  of  rotation. 
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3.0  COMPUTER  PROGRAM 

All  of  the  calculations  necessary  to  predict  the  low  fre- 
quency acoustic  radiation  from  sheet  cavitation  volume  fluctua- 
tions are  done  by  a series  of  interconnected  computer  programs. 

Major  portions  of  the  current  program  were  adapted  from  3BN  pro- 
peller design  and  analysis  programs  developed  under  other  contracts. 

The  program  was  broken  down  into  distinct  modules  to  simplify 
programming  and  debugging.  This  modular  approach  also  allows 
portions  of  the  program  to  be  modified  as  improvements  to  the 
theory  are  made,  without  affecting  other  portions  of  the  program. 

3.1  Computer  Program  Organization 

As  discussed  in  Ref.  1,  the  estimation  of  low  frequency 
acoustic  radiation  from  cavitation  volume  fluctuations  requires 
four  major  steps: 

1)  Select  wakefield; 

2)  Design  prop  to  suit  circumferentially  averaged  wake 
velocities  so  that  proper  thrust  or  horsepower  is  obtained; 

3)  Analyze  prop  in  non-uniform  wakefield  and  determine  the 
volume  of  cavitation  for  each  position  of  the  blade  in 
the  wake; 

4)  Resolve  the  cavitation  volume  vs  time  history  into  a 
Fourier  series  and  compute  the  farfield  radiation  at 
blade  rate  and  harmonics  of  blade  rate. 

The  computer  program  has  been  broken  into  modules  along  these 
same  lines,  as  shown  in  Fig.  3.1.  Each  of  the  major  sections  of 
the  program  are  run  independently,  using  data  files  to  effect  the 
transfer  of  information  between  sections. 
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3.2  Example  Computer  Run 

The  computer  output  showing  the  design  and  analysis  of  the 
baseline  propeller  from  the  sensitivity  analysis  (Section  6.2)  is 
shown  on  the  following  pages.  Fig.  3.2. 
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EFFICIENCY  0.5409  <1-W<V>> 

P'D  R>  0.6335  LAMBDA <MERN>  0 
KT  0.1664  NRCfi  66?fi=0.8  SEC 

BURR ILL  PRRRMETERS:  SIGMfi<.7R> 

WE I GHT  OF  BLADES <LBS> : 1 1 3959 . 0 

REYNOLDS  NUMBER  MULT  I PL I ER=  1 . 0 0 0 0 0 


RPM 

NO.  BLRDES 
EFFICIENCY 
P-D  C.?R> 
KT 


THRUST <LBS>  551343. 
HORSEPOWER  30000. 

SHRFT  IMMER.  50.00 

EXP.RRER  RRTIO  0.600 
< 1 -W  <V>  > 0.607 

LAMBDA  <MEAN>  0.1315 
NRCfi  66 » fl=  0 . 8 SECTIONS 
SIGMA  <.7R>  0.583 


CTS I 1. 
CPS I 1. 
CTS  1 . 
CPS  1 . 
CTHP  6. 
J CMEftN) 
KQ 

TALK.  7R> 


1 . 0843E+00 
1 . 04443+00 
1 . 0739E+00 
1.1 884E+0 0 
6 . 4384E-0 1 
> 0 . 38 1 6 

0. 0184 
R>  0.154 


.V 

TAN  B 

TAN  BI 

GS 

DCT 

DCP 

CLIFT 

CDRAG 

0.30  0 

0. 3448 

0. 7333 

0. 000E+00 

-0. 0005 

0.  0006 

0. 0000 

0. 0070 

0 . 3 0 0 

0. 3669 

0.5367 

3. 343E-03 

0.6002 

0.4917 

0. 4602 

0. 0070 

0 . 4 0 0 

0. 3378 

0.4315 

3 . 656E- 03 

0 . 9638 

0.8573 

0.3555 

0.  0070 

0 . 5 0 0 

0.  3044 

0.3571 

3. 315E-03 

1.3156 

1 . 3673 

0.3813 

0.  007  0 

0 . 6 0 0 

0. 1889 

0.3134 

3.834E-03 

1.6107 

1 . t»y  i j *5* 

0.3345 

0. 0070 

0.  700 

0.  1777 

0.3814 

3 . 796E-03 

1.8801 

3.1300 

0 . 1 869 

0. 0070 

0 . 8 0 0 

0 • 1 668 

0.  3550 

3. 653E-03 

3. 0510 

3.4800 

0. 1591 

0.  007  0 

0 . 9 0 0 

0.  1583 

0.3343 

3. 219E-03 

1 . 9387 

2.5349 

0. 1353 

0. 0070 

1 . 0 0 0 

0. 1514 

0.3173 

0. 000E+00 

0. 0000 

0.  0000 

0. 0000 

0.  0070 

•:S 

VA-'VS 

VT/VS 

Cx'D 

T.-Ti 

T/'C 

FvD  MC 

SKEW 

0 . 2 0 0 

0.345 

0.  000 

0. 1997 

0 . 0468 

0.2345 

0.461 

10.  78 

0 . 3 0 0 

0.40  0 

1 

o 

■ 

Q 

o 

0.  2190 

0 ■ 041 8 

0. 1910 

0.  496 

-8.54 

0 . 4 0 0 

0.455 

-0. 000 

0 . 2364 

0. 0369 

0.  1560 

0.  530 

-6.52 

0.  500 

0.511 

-0. 000 

0.2511 

0. 0319 

0. 1269 

0.  561 

-4 . 78 

0 . 6 0 0 

0. 567 

-0.  000 

0. 2619 

0. 0366 

0.1017 

O.  591 

-3.91 

0 . 7 0 0 

0.  632 

-0.  000 

0. 2663 

0. 0311 

0. 0791 

O.  619 

-0.97 

0 . 8 0 0 

0 .667 

-0.000 

0.3595 

0. 0138 

0. 0530 

O.  641 

1.21 

0 .900 

U ■ 1 2 

-0. 000 

0.3373 

0. 0075 

0.0339 

O.  662 

3.50 

1 . 0 0 0 

0.  757 

0.  000 

0 . 08 1 8 

0. 0033 

0. 0269 

O.  683 

5.51 

v 

/*• 

CLF/CL 

F.'C  2D 

KC 

F/C 

3D  AL  <NT> 

SIGMA 

SIGMA 

0 . 3 0 0 
0 . 3 0 0 
0 . 4 0 0 
0 . 5 0 0 
0 . 6 0 0 
0 . 7 0 0 
0 . 8 0 0 
0 . 9 0 0 
1 . 0 0 0 


1 . 0 0 0 
1 . 0 0 0 
1 . 0 0 0 
1.00  0 
1 . 0 0 0 
1 . 0 0 0 
1 . 0 0 0 
1 . 0 0 0 
1 . 0 0 0 


0. 03969 
U • 03638 
0. 03716 
0.  03091 
0. 01631 
0. 01331 
0. 01107 
0. 00934 
0. 00739 


1 .6083 
1.3136 
1. 0169 
1 . 0004 
1. 0740 
1.1741 
1 . 3593 
1.8133 
3.3671 


0. 04774 
0. 04774 
0. 03761 
0. 03093 
0 . 0 1 753 
0.01 563 
0. 01505 
0 . u 1 6 7 6 
0. 01676 


0. 0000 
1 . 0050 
0.5813 
0. 4403 
0.3687 
0.3389 
0.3168 
0.3533 
0. 0000 


6 . 86 1 1 
3. 9988 
1.6516 
1 . 0353 
0. 7043 
0.  5068 
0 . 38 03 
0.3944 
0. 3335 


0. 0997 
0. 3770 
0.3953 
0.4999 
0.5798 
U. 6337 
0 . 6 085 
0.5530 
0.3977 


IF  CHANGES?  ENTER  I, J? VALUE.  TO  RUN  ENTER  9?  9?  9 66 

INPUT  PROP  GEOMETRY  DATA  FILE  NAME  SRPG 1 
INPUT  DESCRIPTIVE  HERDER?  UP  TO  70  CHARACTERS 
i-R  BASELINE  PROP 

IF  CHANGES?  ENTER  I ?J? VALUE.  TO  RUN  ENTER  9? 9? 9 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
PROM  COPY  PURNISHED  TO  DOC - 

■^OKl  T^?CxSfc<W\  CfWftNL.  "TC>  TfefcS^A  G>C'OC\ ' STC-Titr/ 

Cf\\J  \TPrV\OM  fcV4P\LV£S\£- 

INPUT  PROP  GEOMETRY  BATA  FILE  MR ME  ShPGI 
WAKE  FILE  USED  FOR  PROP  DESIGN  = TRNWK4 


INPUT  TRANSLATED  WAKE  DATA  FILENAME:  TRMWK4 

INPUT  CHANGES  IN  ANALYSIS  CONDITIONS 
INPUT  CHANGE  CODE?  NEW  VALUE:  360 >36 


ENTER  OUTPUT  FILE  HEADER  CUP  TO  65  CHAR> 
A BASELINE  PROPEL  LFR 


ENTER  CAVITATION  ANALYSIS  DATA  FILE  NAME:  SACV1 


THETA= 

0. 

0 

TBLAD=  178619. 

QBLAD=  512883. 

THETA= 

10. 

0 

TBLAD=  172763. 

QELRD=  508368. 

THETA= 

£0. 

0 

TBLAB=  159038. 

QBLAD=  486277. 

THETA= 

30. 

0 

TBLAB=  142613. 

QBLAD=  453980. 

THETA= 

40. 

0 

TBLAD=  128163. 

QBLAD=  422101. 

THET  A= 

50. 

0 

TBLAD=  116900. 

QBLAD=  395579. 

THETA= 

6 0 . 

0 

TBLAD=  108024. 

QBLAD=  373739. 

THETA= 

70. 

0 

TBLAD=  10653. 

BLAD=  354975. 

THETA= 

8 0 . 

0 

TELAD= 

94598. 

QBLAD=  333998. 

THETA= 

9 0 . 

0 

TBLRB= 

89827. 

QBLAD=  325956. 

THET  A= 

1 0 0 . 

0 

TBLAD= 

86504. 

QBLAD=  316367. 

THETA= 

110. 

0 

TBLAD= 

85064. 

QBLAD=  311409. 

THETA= 

1 2 0 . 

0 

TELAD= 

8671 8 . 

QBLAD=  31 4241.  *2— 

THETA= 

130. 

0 

TBLAD= 

93141. 

QBLAD=329054. 

THETA= 

140. 

0 

TBLAD*  104961. 

QELAD=  356807. 

THETA= 

150. 

0 

TBLAD=  120210. 

QELAB=  392274. 

THETA= 

1 6 0 . 

0 

TBLAD=  131518. 

QBLAD=  417763. 

THETA= 

170. 

0 

TBLRD=  132500. 

QBLRD=  419987. 

THETA= 

180. 

0 

TBLAD=  121492. 

QBLAD=  395118. 

THETA= 

190. 

0 

TBLAD=  103722. 

QBLRD=  352560. 

THETA= 

£ 0 0 . 

0 

TBLAD= 

86174. 

QBLAD=  307433. 

THETA= 

£10. 

0 

TBLAD= 

72933. 

QBLAD=  271159. 

THETA= 

££0. 

0 

TBLAD= 

65370. 

QBLAD=  249184. 

THETA= 

£30. 

0 

TBLAD= 

62450. 

QBLAD=  £39328. 

THETA= 

£4  0. 

0 

TBLAD= 

62622 . 

QBLAD=  239150. 

THETA= 

£50. 

0 

TBLAD= 

64569. 

QBLAD=  £43504. 

THETA= 

£6  0 . 

0 

TBLAD= 

67725. 

QBLAD=  £51209. 

THETA= 

£70. 

0 

TBLAD= 

72320. 

QBLAD=  £62712. 

THETA= 

280. 

0 

TBLAD= 

78700. 

QBLAD=  £78677. 

THETA= 

£90. 

0 

TBLAD= 

87237. 

QELRB=  £99760. 

THETA= 

3 0 0 . 

0 

TBLAD= 

98480. 

QBLAD=  327007. 

THETA= 

3 1 0 . 

0 

T BLAD=  113376. 

QBLAD=  362247 . 

TH£TA= 

3£  0 . 

0 

TBLAD=  131284. 

QBLAD=  403067. 

THETA= 

330. 

0 

TBLAD=  150149. 

QBLAD=  444376. 

THET  A= 

340. 

0 

TBLAD=  165973. 

QBLAD=  473299. 

THE  TA  = 

350. 

0 

T BLAD=  176122. 

QBLAD=  501971. 

Outpcts*  Tneo'sr 

T Tc5v?G>oe  Pcte. 
"BuAtie;  fcr  cacm 

& 
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LOW  FREQUENCY  LI  HE  CALCULATION 

♦♦♦  ♦♦♦♦♦♦♦♦♦  ♦♦♦♦  ♦♦♦♦♦♦♦♦♦♦♦ 


CAVITATION  ANALYSIS  DATA  FILE: 
SA  BASELINE  PROPELLER 

WAKE  DATA  FILE:  TANWK4 

PROP  GEOMETRY  DATA  FILE:  SAPG1 


SAC  VI 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
FROM  copy  furbished  to  DOC  " 


VSHIP  <KTS>  16.00 

D I AM  (FT) 

30.  00  SHAFT  IMME 

RPM 

■36.  00 

NO.  BLADE 

S 5 EXP.  AREA 

THRUST 

543405. 

KT 

0.1641  SIGMA  CVS) 

DHP 

£94£4. 

1 0*KQ 

0.  1808  SI GMA  <N> 

VALUES 

FOR  INDEX  BLADE: 

THETA 

CAVITY  AREA 

AC''AB 

CAVITY  VOLUME 

DEG 

<FT»-*£> 

<FT**3> 

0 . 0 0 

8.59 

0.  105 

0.54685 

1 0.  0 0 

9.  £8 

0.114 

0.  50992  C ~v2fW 

£0.  00 

£.37 

0.029 

0 . 06488  C 

30.  00 

0.  36 

0.  004 

0. 00221 J 

40.  00 

0.  00 

0.  000 

0. 00000 

5 0 . 0 0 

0.  00 

0.  000 

0. 00000 

60.  00 

0 . 0 0 

0.  000 

0. 00000 

70.  00 

0.  00 

0.  000 

0. 00000 

SO.  00 

0.  00 

0.  000 

0. 00000 

90.  00 

0.  00 

0.  000 

0. 00000 

1 0 0 . 0 0 

0.  00 

0 .000 

0. ooooo 

1 1 0.  00 

0.  00 

0.  000 

0. ooooo 

1 £ 0 . 0 0 

0.  0 0 

0.  000 

o. ooooo 

1 3 0 . 0 0 

0.  00 

0.  000 

o. ooooo 

140.  00 

0 . 0 0 

0.  000 

o. ooooo 

150.  00 

0.  00 

0.  000 

o. ooooo 

160. 0 0 

0.  00 

0.  000 

o. ooooo 

170.  0 0 

0 . 0 0 

0.  000 

o. ooooo 

190.  00 

0.  00 

0.  000 

o. ooooo 

190.0  0 

0.  00 

0.  000 

o. ooooo 

£ 0 0 . 0 0 

0.  00 

0.  000 

o. ooooo 

£10. 00 

0.  00 

0.  000 

o. ooooo 

££  0 . 0 0 

0.  00 

0.  000 

0. ooooo 

£3 0 . 00 

0.  00 

0.  000 

o. ooooo 

£40. 00 

0.30 

0.  004 

0.  00239 

250.  00 

0.  31 

0 . 0 04 

0.00239  ( 

260.  00 

0.13 

0.  002 

0.00082  v 

£70.  00 

0.  14 

0.  002 

0.00082  1 

£80.  00 

0 . 0 0 

0.  000 

0.00  0 0 0 

£90.  00 

0.  00 

0.  000 

o. ooooo 

300.  00 

0.  00 

0.  000 

0 . 0 0 0 0 0 

310.  0 0 

0.  00 

0.  000 

o. ooooo 

320. 00 

0.  00 

0.  000 

o. ooooo 

330.  00 

0.  79 

0.  010 

0.  01052  **} 

34  0.  00 

3.79 

0.  047 

0.  14180  s 

35  0 . 0 0 

7.59 

0.  093 

0.41666  J 

AVERAGE 

CAVITY  VOLUME 

<FT*4-3>  : 

0.23654  (FC*2. 

50.  00 
0.60  0 
7.271 
£.872 


CAV  \1ATvOt4 


(T  fTT\ 
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SACVt 

<B('asuKie.  PP<b9€u.eK 


=1 


LINE  LEVEL  CALCULATIONS:  HARMONICS  OF  BLADE  RATE 
ACOUSTIC  PRESSURES  IN  DB  RE  1 MICROPASCAL  AT  1 YARD 


N 

FREQ 

A <N  > 

MONOPOLE 

LLOYD 

LRESL 

MONOPOLE 

•%FT  ♦♦3.> 

SPL  <DB> 

MIRROR  <BB> 

<DB> 

S' 30  FT  <DB> 

1 

7.  17 

0. £99664 

180.7 

-20.4 

160.3 

185.2 

2 

14.33 

0. 058664 

1 78 . 6 

-14.5 

164.2 

183.  1 

3 

£1.50 

0. 033294 

180.7 

-11.0 

169.7 

185.2 

4 

£3.67 

0.  01 54  0 0 

179.  0 

—8. 6 

170.5 

183.5 

5 

35 . 33 

0. 003401 

169.8 

-6.7 

1 63 . 1 

174.2 

6 

43.  00 

0. 00£303 

169.6 

—5. 3 

164.3 

174.  0 

1 

50.  17 

0. 000448 

158.  0 

-4.  1 

153.9 

162.5 

o 

57.33 

0. 000594 

162.8 

-3.  1 

159.7 

167.3 

9 

64.50 

0. 000846 

167.9 

-2.3 

165. 6 

172.4 

10 

71.67 

0. 001235 

173.  0 

— 1.6 

171.4 

177.5 

END  OF  EXECUTION 
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4.0  EXAMPLE  CALCULATIONS  USING  MERCHANT  SHIP  PROPELLER  PROGRAM 

Several  examples  of  calculations  done  using  the  current  com- 
puter programs  are  presented  in  this  section  to  illustrate  the 
accuracy  and  the  shortcomings  of  the  program. 

4.1  Calculated  and  Observed  Cavitation  Patterns 

4.1.1  Single  screw  merchant  ship 

The  first  example  is  taken  from  ref.  14  and  deals  with  the 
propeller  for  a medium  size  single  screw  cargo  vessel,  with  the 
following  particulars: 

Ship  speed  = 15.8  knots 
Prop  diameter  = 15.9  feet 
RPM  = 155.5 

Shaft  horsepower  = 9500 
No.  of  blades  = 4 

Expanded  Area  Ratio  = 0.6  (ratio  of  blade 

area  to  propeller 
disc  area) 

The  cavitation  properties  of  this  propeller  were  determined  by 
running  a model  in  a cavitation  tunnel  in  a simulated  ship  wake- 
field.  The  results  of  the  cavitation  tunnel  tests  are  shown  in 
Pig.  4.1,  along  with  the  predictions  computed  in  ref.  14  and  the 
BBN  predictions  using  the  current  computer  programs. 

Van  Cassanen's  predicted  cavitation  patterns  come  much  closer 
to  matching  the  observed  cavitation  patterns.  Note,  however,  that 
van  Oossanen's  predictions  are  for  the  model  propeller  while  the 
BBN  predictions  are  for  the  full  scale  propeller.  Reference  14 
points  out  that  the  large  Reynolds  number  difference  between  model 
and  full  scale  propellers  can  have  considerable  impact  on  the  cav- 
itation patterns  experienced  under  otherwise  identical  flow  con- 
ditions. The  theory  developed  by  van  Oossanen  in  ref.  14  goes  to 
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great  lengths  to  account  for  the  effect  of  Reynolds  number  on 
cavitation,  which  causes  the  theory  to  be  much  more  complicated 
than  the  one  developed  by  BBN.  However,  both  theory  and  experi- 
ment show  that  the  methods  used  in  the  BBN  cavitation  prediction 
scheme  are  approximately  correct  for  very  high  Reynolds  numbers; 
i.e.,  full  scale  ship  propellers.  The  fact  that  the  BBN  predic- 
tion shows  sheet  cavitation  on  the  pressure  side  while  the  model 
tests  do  not  is  in  line  with  full  scale  experience:  this  form  of 
cavitation  has  been  observed  on  full  scale  propellers  even  though 
none  was  present  during  model  testing.  The  extent  of  sheet  cav- 
itation on  the  suction  side  of  the  blade  may  also  change  with 
Reynolds  number. 

The  major  shortcoming  of  the  BBN  cavitation  prediction  scheme 
is  the  lack  of  a way  to  handle  cavitation  in  the  tip  vortex  which 
comes  off  of  the  tip  of  the  blade  (similar  to  an  aircraft  tip  vor- 
tex). There  are  as  yet  no  simple  methods  for  predicting  the  volume 
and  extent  of  tip  vortex  cavities,  but  it  appears  that  the  cavita- 
tion volume  associated  with  the  tip  vortex  on  merchant  ship  pro- 
pellers will  be  small  compared  to  the  volume  resulting  from  sheet 
cavitation.  Hence  the  tip  vortex  cavitation  volume  may  influence 
the  generation  of  noise  at  the  higher  harmonics  of  blade  rate, 
but  probably  won't  affect  the  acoustic  radiation  at  the  lower 
harmonics,  which  are  of  primary  interest  in  the  present  work. 

4.1.2  Twin  screw  merchant  ship 

The  second  example  is  also  taken  from  van  Oossanen's  work 
(ref.  14)  and  concerns  the  propeller  of  a fast  twin  screw  vessel 
with  the  following  particulars: 

Design  ship  speed  = 30  knots 
Prop  diameter  = 11. 98  feet 
Design  RPM  = 230 
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Design  horsepower  = 15500 
No.  of  blades  = 5 
Expanded  Area  Ratio  = 0.8 

The  observations  and  predictions  shown  in  Fig.  4.2  are  for 
a ship  speed  of  24  knots,  with  the  prop  turning  at  176  RPM.  It 
is  seen  that  the  correlation  between  full  scale  observations  and 
the  BBN  prediction  is  satisfactory,  with  the  exception  of  the  tip 
vortex  cavity. 

It  is  interesting  to  note  that  a model  test  of  the  same  pro- 
peller under  identical  flow  conditions  (except  for  the  smaller 
Reynolds  number  associated  with  the  model)  did  not  show  any  evi- 
dence of  sheet  cavitation  on  the  pressure  side  of  the  blade,  as 
is  shown  in  Fig.  4.2  (283°  and  312°  positions).  This  reinforces 
the  discussion  in  Section  4.1.1  concerning  Reynolds  number  effects 
on  cavitation  patterns. 

4.2  Shortcomings  of  Prediction  Process 

It  appears  that  the  major  problem  with  the  BBN  cavity  predic- 
tion process  is  the  neglect  of  the  tip  vortex  cavity.  As  noted 
previously,  neglecting  this  phenomenon  will  probably  only  affect 
the  radiated  noise  at  the  higher  harmonics  of  blade  rate,  so  that 
it  would  appear  to  be  safe  to  neglect  this  in  the  current  work 
which  is  mainly  concerned  with  correctly  modelling  the  first  few 
harmonics  of  blade  rate.  The  effect  of  the  tip  vortex  cavity  will 
be  included  as  soon  as  time  and  funds  permit. 

A brief  note  on  the  whole  concept  of  mathematically  modelling 
propeller  cavitation  is  in  order  here.  Although  it  is  clear  that 
the  prediction  process  is  subject  to  continual  improvement  and 
refinement,  there  is  a danger  that  the  associated  calculations 
can  get  so  complicated  and  expensive,  and  require  so  much  input 
data,  that  they  are  of  little  use  in  solving  the  problem  at  hand 
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(predicting  ocean  ambient  noise).  We  believe  that  the  present 
model  is  probably  adequate  for  predicting  radiated  noise  at  the 
first  few  harmonics  of  blade  rate  but  will  require  some  improve- 
ments, such  as  including  the  tip  vortex  cavity,  before  it  is  use- 
ful for  predicting  the  cavitation  noise  at  the  higher  harmonics 
of  blade  rate. 

Further  comparisons  (beyond  the  two  presented  here)  of  pre- 
dicted cavitation  patterns  with  those  measured  on  full  scale  pro- 
pellers should  be  conducted  before  exercising  this  program  to  pre- 
dict the  ensemble  characteristics  of  ocean  ambient  noise  due  to 
merchant  ship  propeller  cavitation.  These  comparisons  are  neces- 
sary in  order  to  establish  whether  or  not  the  current  program 
provides  unbiased  estimates  of  propeller  cavitation  volume. 
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5.0  REVIEW  OF  LINE  LEVEL  PREDICTION  PROCESS 

5.1  Radiation  from  Cavitation  Volume  Fluctuations 

For  the  purposes  of  calculating  acoustic  radiation,  the  pre- 
sent work  assumes  that  the  sheet  cavitation  volume  fluctuations 
may  be  adequately  represented  by  an  equivalent  monopole  at  an 
appropriate  depth  below  the  surface  (currently  assumed  to  be  at 
the  propeller  shaft  centerline).  Then  the  pressure  radiated  from 
this  simple  volume  source  is: 

<p2>  = <Q2>  (5.1) 

where  <p2>  = mean  square  pressure 
Q = volume  velocity 
f = frequency  (Hz) 
r = radius  from  acoustic  enter 
p = density  of  fluid 

If  we  represent  the  cavitation  volume  history  for  the  whole  pro- 
peller as  a Fourier  series 

00  IN  tot 

Q(t)  = l A e (5.2) 

N = 0 W 

where  Q(t)  = cavitation  volume  on  prop 
= complex  coefficient 

N = harmonic  number 


go  = blade  rate  frequency  (rad/sec)  = RF^Pr0P  • 2ttB 
B = number  of  propeller  blades 
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and 


<Q2  > = i .X 


:Q2>  " f N^O  N2-2IAnI2 


(5.4) 


Therefore  the  radiated  pressure  at  the  N-th  harmonic  of  blade 
rate  is 

r2,„2  | A. 


<P2>  = ^ 


(5.5) 


The  frequency  f^  of  the  N-th  harmonic  of  blade  rate  is 

r = Nw 
N 2 it 


(5.6) 


So  that  Eq.  5.5  reduces  to 


<P2>  = 2 it2  fH2  |AN| 2 


N 

• N ' I AN 


TT2f„4|A„|2p2 


(5.7) 


2 r2 


The  problem  of  predicting  the  cavitation  line  spectrum  at  blade 
rate  and  harmonics  thus  consists  of  finding  the  coefficients  A^ 
in  Fourier  series  expansion  of  the  propeller  cavitation  volume 
history . 

The  coefficients  A^  may  be  determined  from  the  volume  history 
of  a single  propeller  blade,  as  shown  below.  Representating  the 
cavitation  volume  on  the  first  (index)  blade  as  a Fourier  series 
expansion  in  rotation  angle  6(0<6<2it) 

V (0)  = l C ein6  (5.8) 

1 n=0  n 

where  C = complex  coefficient 
n 

and  assuming  all  of  the  blades  have  identical  volume  histories, 
the  volume  on  the  bth  blade  will  be 

Vb(6)  . V,  (0  + 

and  the  volume  history  for  the  whole  propeller  will  be 


l 

J 
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n=0 
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b = l 


in(0+^) 


C e 
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oo 
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n=0 


C e 
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2 7Tb 

B 


(5.9) 


The  first  summation  above  is  the  volume  history  of  the  index  blade 
[V  (0)]  and  the  second  summation  has  nonzero  values  only  for  values 
of  n equal  to  integer  multiples  of  B,  i.e.. 


Y in  B B,  n=k*B 

I e = 

b=l  0,  n/k*B 


(5.10) 


Thus  the  cavitation  volume  history  on  a cavitating  propeller  with 
B blades,  in  terms  of  the  Fourier  series  expansion  coefficients 
for  the  volume  history  on  the  index  blade,  Cn,  is 


V (9) 
P 


B 


OO 


l 

k=0 


i (Bk0 ) 

B *K  e 


(5.11) 


Since  rotation  angle  0 and  time  t are  directly  related  by  the  ro- 
tational speed  of  the  propeller,  with  the  period  of  rotation  T 
corresponding  to  2tt  radians  of  rotation,  we  may  write 


yo 

Comparing  equations  (5.2) 
so  that 


oo 


= B 

and 


l C 

k=0  B*k 
(5.12), 


. 2irBk 

m 


we  see  that  Q(t) 


(5.12) 


V (t) 
P 
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I 

N=0 


V 


iNu)t 


= B 


l 

k=o 


2irBk 


'B*k 


(5.13) 


Changing  the  index  of  summation  k to  N on  the  right  hand  side, 

and  noting  that  2t:B  = to  (blade  rate  frequency),  we  obtain 

T 


I A eltoNt  = b l C ela)Nt 
N=0  N N=0  B*N 


so 

or 


A = B*Cn  .. 
N B *N 


'AN'  B ’ I CB -N I 


(5.14) 


Thus  to  get  the  radiated  acoustic  pressure  at  the  Nth  harmon- 
ic of  blade  rate,  we  must  compute  the  (B*N)th  harmonic  of  the 
index  blade  cavitation  volume  history. 

This  immediately  causes  problems  in  predicting  the  radiated 
noise  at  higher  harmonics  of  blade  rate,  since  high  harmonics  of 
the  index  blade  volume  history  are  required.  Although  the  magni- 
tude of  these  higher  harmonics  will  be  small,  equation  (5.7)  shows 
that 

<p!>~VIanI2 


so  that  the  f^1*  term  will  cause  the  higher  harmonics  (up  to  about 
7th  and  possibly  beyond)  to  radiate  significant  power. 


5.3  Method  of  Calculating  A 


N 


The  current  computer  program  calculates  the  cavitation  volume 
on  the  index  blade  at  evenly  spaced  values  of  6.  The  choice  of 


• » 
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the  spacing  is  a compromise  between  two  opposing  considerations: 
computer  costs  and  accuracy.  For  this  sensitivity  analysis,  most 
of  the  computations  have  been  done  using  10°  Increments  of  0, 
yielding  36  data  points  (360°/10°  = 36).  If  we  were  to  do  a har- 
monic analysis  on  only  these  36  original  points,  then  we  could 
calculate  only  the  first  18  harmonics  at  most  (Nyqulst  frequency). 
For  a 5 bladed  propeller,  Eq.  (5.14)  tells  us  that  we  could  then 
compute  the  radiated  noise  only  up  to  the  third  harmonic 
( | A 3 [ = 5 |Cj5|)>  Furthermore,  since  there  may  be  frequency 
content  beyond  twice  the  folding  frequency,  there  may  be  aliasing 
even  below  the  Nyquist  frequency. 

To  attempt  to  alleviate  some  of  the  aliasing  problems,  the 
computer  program  is  set  up  to  pass  a "smooth"  curve  (continuous 
first  derivative)  through  the  calculated  points  and  interpolate 
additional  points  (for  a total  of  240  points)  before  performing 
the  harmonic  analysis.  Obviously,  this  procedure  is  only  as  valid 
as  the  assumption  that  the  growth  and  collapse  of  the  sheet  cav- 
itation bubble  is  a "smooth"  process.  This  points  is  discussed 
further  in  Section  6.4. 
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6.0  SENSITIVITY  ANALYSIS 

6.1  Purpose 

The  purpose  of  this  sensitivity  analysis  is  to  determine  the 
influence  of  analysis  parameters,  propeller  characteristics,  and 
wakefield  characteristics  on  the  prediction  of  cavitation  line 
spectra.  This  will  provide  guidance  in  improving  the  prediction 
procedure  and  suggest  which  factors  are  important  in  influencing 
the  cavitation  noise  from  merchant  ships.  It  is  these  factors 
which  must  be  determined  for  the  world's  merchant  ship  fleet  in 
order  to  estimate  the  ensemble  characteristics  of  the  ocean  ambient 
noise  due  to  merchant  ship  propeller  cavitation. 

6 . t Baseline  Propeller  and  Wakefield 

The  baseline  propeller  for  this  sensitivity  analysis  is 
similar  to  a propeller  which  might  be  found  on  a 300,000  deadweight 
ton  tanker.  The  wakefield  for  such  a ship  was  estimated  from  the 
data  presented  in  Ref.  17,  and  is  shown  in  Figures  6.1  and  6.2. 

The  baseline  propeller  was  designed  to  match  this  wake  and  has 
the  following  characteristics: 

Diameter  = 30.0  ft 
Design  RPM  = 86 

Design  ship  speed  - 16.0  knots 

Horsepower  at  design  speed  = 30,000  SHP 

Number  of  blades  = 5 

Expanded  Area  Ratio  = 0.6 

Percentage  of  lift  by  camber  = 100 % 

Immersion  to  shaft  g = 50.0  ft 

Optiirum  pitch  distribution 
(maximum  efficiency) 

B-Series  skew  distribution 

Blade  thickness  based  on  strength  criteria 


"Example  Wfvte  vas.  CU=.©3  hpimkcr.  csf  OncnouJn 
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6.3  Number  of  Analysis  Points 

A necessary  check  on  any  sort  of  numerical  calculation  is  a 
convergence  check.  The  present  computer  program  appears  to  have 
no  problems  in  converging  to  a cavity  volume  for  a given  blade 
position  in  the  wake.  The  real  question  is  how  many  analysis 
points  (how  many  0 values)  are  required  to  get  reliable  results 
for  the  higher  harmonics.  In  order  to  get  a handle  on  this  prob- 
lem, the  baseline  propeller  was  analyzed  in  the  baseline  wake  for 
3 different  0 spacings:  10°,  5°,  and  3°,  yielding  36,  72,  and  120 
points,  respectively.  Note  that  the  computer  program  still  inter- 
polates additional  points  before  performing  the  harmonic  analysis, 
as  described  in  Section  5-3.  The  predicted  monopole  source  levels 
at  the  shaft  centerline,  in  dB//(l  pPaxYD)2,  are  shown  in  Fig.  6.3 
for  these  3 analysis  spacings.  It  is  obvious  that  the  analysis 
spacing  becomes  increasingly  important  as  the  harmonic  number  in- 
creases. The  reason  for  the  large  discrepancies  at  the  higher 
harmonics  of  blade  rate  may  be  seen  from  Fig.  6.4,  which  show  the 
growth  and  collapse  of  the  sheet  cavitation  volume  as  the  blade 
passes  through  the  low  velocity  region  behind  the  ship  in  the 
vertical  upright  position.  The  curves  shown  are  the  "smooth" 
curves  which  the  computer  program  passes  through  the  calculated 
points  in  order  to  interpolate  further  points  for  the  harmonic 
analysis,  as  described  in  Section  5.3.  It  is  obvious  that  rela- 
tively small  spacing  of  points  is  necessary  in  order  to  define 
all  of  the  lumps  and  bumps  in  the  curve,  which  greatly  affect  the 
higher  harmonics.  The  major  reason  for  not  using  a fine  spacing 
all  of  the  time  is  the  cost  of  the  computer  calculations.  The 
approximate  cost  of  the  analysis  for  a given  prop  in  a given  wake 
is  shown  below,  as  a function  of  analysis  point  spacing: 


‘SEMsmMiTY  anAlvsis.  - EFFecrr  GP  Am^LY'S.S 
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Analysis  Spacing  Number  of  Points  Approx  Cost 


10° 

36 

$14 

5° 

72 

$27 

3° 

120 

$45 

While  these  costs  are  reasonable,  it  is  obvious  that  the  total 
bill  can  become  substantial  if  many  propellers  are  analyzed  with 
a very  small  spacing.  Since  it  is  not  yet  clear  that  the  higher 
harmonics  may  be  accurately  calculated  without  some  consideration 
of  bubble  dynamics  and  the  tip  vortex  cavity,  the  sensitivity 
analysis  presented  here  was  done  using  10°  analysis  spacing 
for  the  most  part.  Only  the  first,  second,  and  third  harmonics 
of  blade  rate  will  be  reported  for  the  majority  of  the  cases. 

6.4  Generation  of  Higher  Harmonics 

In  deciding  whether  or  not  the  quasi-steady  type  of  cavita- 
tion analysis  used  in  the  present  work  is  really  useful  in  pre- 
dicting the  cavitation  noise  at  higher  harmonics  of  blade  rate, 
it  is  necessary  to  determine  the  cause  of  these  higher  harmonics. 
There  are  three  major  possible  causes;  two  physical  and  one  re- 
lated to  numerical  analysis  problems.  These  are  listed  below. 

o The  response  (change  in  loading)  of  the  propeller  blade  to 
the  higher  harmonics  of  the  wakefield. 

o The  non-linearities  and  threshold  phenomena  present  in 
the  cavitation  process  (non-linear  response  of  a system 
to  a sinusoidal  input  will  generate  harmonics). 

o The  discretization  of  the  blade  when  calculating  the  cav- 
itation volume  present. 
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In  order  to  see  which  of  these  phenomenae  are  important,  the  base- 
line propeller  was  analyzed  in  a perfectly  sinusoidal  wakefield 
(one-'  per  revolution  sinusoidal  variation  in  inflow  velocity). 

The  resulting  predicted  source  levels  at  the  various  harmonics  of 
blade  rate  are  shown  in  Pig.  6.5  It  is  clear  that  response  of  the 
blade  to  higher  harmonics  of  the  inflow  perturbations  is  not  nec- 
essary to  generate  high  harmonics,  and  implies  that  the  response 
of  the  blade  need  only  be  calculated  for  the  first  few  harmonics 
of  the  inflow  perturbations.  This  is  very  encouraging,  since  the 
quasi-steady  lifting  line  theory  used  in  the  present  computer  pro- 
gram is  perfectly  adequate  for  performing  these  calculations.  If 
the  response  of  the  blade  to  the  higher  harmonics  at  the  inflow 
perturbations  were  required,  a far  more  complicated  unsteady  lift- 
ing surface  theory  would  be  needed  to  correctly  calculate  the  re- 
sponse of  the  blade  to  inflow  perturbations. 

A look  at  the  actual  volume  vs  0 history  for  this  propeller 
(Pig.  6.6)  is  quite  revealing.  The  jumps  in  the  volume  at  0 
values  of  approximately  310°  and  30°  are  indeed  due  to  numerical 
discretization  errors.  These  arise  because  the  cavity  cross  sec- 
tional areas  which  are  integrated  in  the  radial  direction  to  cal- 
culate cavity  volume  are  calculated  only  at  a finite  number  of 
radial  positions  on  the  blade.  The  sudden  jump  in  the  calculated 
volume  is  due  to  a sudden  calculated  jump  in  the  radial  extent  of 
cavitation.  In  reality,  though,  the  radial  extent  of  the  sheet 
cavitation  bubble  will  tend  to  vary  smoothly.  Hence  this  sudden 
jump  is  mostly  a problem  associated  with  the  numerical  procedure 
used.  These  jumps  or  discontinuities  in  the  volume  history  will 
contribute  to  the  higher  harmonics. 

Note,  however,  that  even  if  the  discontinuities  mentioned 
above  were  not  present  in  Pig.  6.6,  there  would  still  be  signifi- 
cant higher  harmonic  content.  This  may  be  illustrated  as  follows. 
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The  calculated  volume  history  may  be  represented  (at  least 
to  first  order)  by  a triangular  shape,  as  shown  in  Pig.  6.7.  To 
get  the  resulting  frequency  content,  we  take  the  Fourier  transform 
of  this  triangular  distribution,  which  is 


sin2  ^ 


F(w)  = 


(6.1) 


A graph  of  amplitude  (F(uO)  vs  frequency  (oj)  for  this  function 
will  obviously  have  nulls  because  of  the  sin2  term.  However,  if 
we  restrict  our  attention  only  to  the  upper  envelope  of  the  func- 
tion and  concern  ourselves  only  with  the  dependence  on  frequency 
(in  order  to  bound  the  problem),  we  may  write 


F(w)  ~ ~ 
, .2 


(6.2) 


Referring  to  the  radiation  relationships  developed  in  Section  5.1, 
we  see  that  for  an  observer  at  a constant  distance  from  a volume 


source 


<p2>  ~ VIV2 


(6.3) 

( from  Eq . (5.7)) 


Noting  that 


fN  ~ w 


|An|  ~ F(w) 

then  substituting  into  Equation  6.3  yields 


FIGURE  6.7 
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~ Constant  (6.4) 

So  the  upper  envelope  of  the  monopole  source  strengths  at  harmon- 
ics of  blade  rate  will  tend  to  be  constant  with  frequency  for  the 
triangular  volume  history  assumed.  Since  real  volume  histories 
will  not  have  the  sharp  corners  of  this  assumed  triangular  shape 
the  higher  harmonics  will  fall  off  eventually. 

The  above  discussion  and  example  demonstrate  that  the  physical 
basis  for  the  higher  harmonics  of  blade  rate  is  the  nonlinear, 
threshold  relationship  between  blade  loading  and  cavity  volume. 

6.5  Ship  Speed  Reduction 

The  effect  of  ship  speed  reduction  on  the  cavitation  source 
strength  is  shown  in  Fig.  6.8  for  the  baseline  propeller  and  wake- 
field.  The  effect  is  certainly  substantial,  amounting  to  a 30  dB 
reduction  in  source  level  for  a 30%  reduction  in  ship  speed.  This 
means  that  the  speed  at  which  ships  are  running  (as  a percentage 
of  design  speed)  must  be  considered  if  ocean  ambient  noise  from 
merchant  shipping  is  to  be  properly  estimated. 

6.6  Blade  Number 

The  effect  of  blade  number  on  monopole  source  level  is  small 
for  "equivalent"  propellers,  as  shown  by  Fig.  6.9.  For  this  sensi- 
tivity analysis  "equivalent"  propellers  have  been  defined  as  hav- 
ing the  same  margin  against  back-bubble  cavitation,  a form  of 
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cavitation  which  generally  causes  severe  erosion  of  propeller 
blades.  This  method  of  defining  equivalent  propellers  is  in  line 
with  standard  naval  architectural  practice. 

Although  the  monopole  source  strength  does  not  vary  signifi- 
cantly with  blade  number,  the  variation  in  blade  number  at  constant 
propeller  RPM  implies  a variation  in  blade  rate  frequency  which 
will  change  the  Lloyd  mirror  and  consequently  the  pressure  propa- 
gated to  long  ranges. 

6.7  Expanded  Area  Ratio 

The  expanded  area  ratio  (EAR)  is  defined  as  the  ratio  of 
the  propeller  blade  area  to  the  propeller  disc  area,  and  hence 
is  just  a non-dimensional  measurement  of  blade  area.  The  choice 
of  EAR  for  a given  propeller  is  a compromise  between  efficiency 
and  cavitation  considerations.  The  EAR  should  be  as  small  as  pos- 
sible to  minimize  the  viscous  drag  on  the  blades  and  hence  maxim- 
ize propeller  efficiency.  However,  if  the  EAR  is  too  small,  cav- 
itation will  be  excessive,  leading  to  excessive  noise  and  vibra- 
tion and  finally,  thrust  breakdown. 

Figure  6.10  indeed  shows  that  the  effect  of  EAR  on  cavitation 
is  quite  large.  Hence  this  variable  must  be  properly  modeled  in 
order  to  predict  the  ocean  ambient  noise  due  to  propeller  cavita- 
tion. 

6.8  Percentage  Lift  by  Camber 

A propeller  is  designed  to  suit  the  mean  (circumferential 
average)  wakefield.  The  design  process  consists  of  calculating 
the  flow  angle  and  lift  coefficient  at  each  radius  so  that  the 
desired  thrust  or  power  absorption  is  obtained.  The  desired  lift 
coefficient  may  be  obtained  by  either  cambering  the  blade  or  pro- 
viding a flat  blade  with  an  angle  of  attack  (see  Fig.  2.3),  or 
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some  linear  combination  of  the  two.  The  highest  propeller  effi- 
ciency will  generally  be  obtained  if  all  of  the  lift  is  generated 
by  camber,  but  in  a non-uniform  wakefield  a highly  cambered  blade 
(100%  lift  by  camber)  may  experience  sheet  cavitation  on  the  pres- 
sure side  of  the  blade,  as  depicted  in  Fig.  4.1.  This  form  of 
cavitation  causes  blade  erosion  and  is  avoided  whenever  possible. 

One  method  of  eliminating  face  cavitation  is  to  design  the  propel- 
ler so  that  some  of  the  mean  lift  is  generated  by  angle  of  attack, 
thus  reducing  the  percentage  of  lift  by  camber. 

Figure  6.11  shows  that  varying  the  percentage  lift  by  camber 
has  a significant  effect  on  cavitation  noise.  Since  the  baseline 
propeller  (100%  lift  by  camber)  is  predicted  to  have  slight  face 
cavitation,  the  prop  with  90%  lift  by  camber  is  probably  closer 
to  the  propeller  that  would  be  found  on  a large  tanker  with  a 
wakefield  similar  to  the  one  used  for  the  sensitivity  analysis. 

6.9  Pitch  Distribution 

The  radial  distribution  of  propeller  blade  pitch  effects 
both  the  efficiency  and  cavitation  performance  of  marine  propellers. 
Current  wisdom  is  to  slightly  reduce  the  pitch  at  the  tip  of  the 
blade  in  order  to  reduce  cavitation.  Figure  6.12  shows  that  this 
approach  does  work,  but  that  the  effect  is  relatively  small. 

6.10  Propeller  Diameter 

I Geometrically  similar  cavitation  on  different  size  pro- 
pellers, the  prop  with  the  larger  diameter  will  have  a greater 
cavitation  source  strength  because  the  volume  of  cavitation  will 
be  proportional  to  (diameter)3.  Since  the  amount  of  cavitation 
present  on  merchant  ship  propellers  at  full  speed  tends  to  be 
constant  fraction  of  blade  area  (see  discussion  in  Sec.  6.7),  it 
is  clear  that  propeller  diameter  is  a major  variable  affecting  the 
cavitation  volume  present  on  a marine  propeller,  and  hence  the 
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acoustic  source  strength. 

6.11  Des i gn  RPM 

The  effect  of  design  RPM  on  "equivalent"  propellers  is  fairly 
small,  as  shown  by  Fig.  6.14. 

6.12  Design  Shaft  Horsepower  (SHP) 

Changing  the  design  SHP  while  holding  a constant  margin 
against  back  bubble  cavitation  results  in  a comparatively  small 
change  in  monopole  source  level.  This  is  shown  in  Fig.  6.15. 

6.13  Propeller  Submergence 

For  a given  propeller  and  wakefield,  decreasing  the  depth  of 
submergence  will  increase  the  amount  of  cavitation  and  thus  the 
monopole  source  strength.  For  the  example  calculated  for  Fig. 
6.16,  decreasing  the  immersion  to  shaft  centerline  from  50  feet 
to  40  feet  caused  the  blade  rate  source  level  to  increase  by  5 dB. 
Note,  however,  that  this  is  for  the  same  propeller  and  wakefield. 
If  the  prop  were  designed  to  operate  at  a shaft  immersion  of  40 
feet  then  the  expanded  area  ratio  would  have  to  be  increased  to 
keep  the  margin  against  back  bubble  cavitation  the  same,  resulting 
in  only  a minor  change  in  radiated  noise.  Hence  the  design  sub- 
mergence is  not  very  important  in  determining  cavitation  noise 
(except  for  the  change  in  Lloyd  mirror),  but  the  operating  sub- 
mergence of  a given  propeller  can  be  important.  This  depends  on 
whether  a ship  is  operating  in  the  loaded  condition  or  in  the 
ballast  condition. 

6.14  Wakefields 

The  baseline  propeller  was  analyzed  in  several  different 
wakefields,  each  of  which  has  the  same  circumferential  average 
velocities,  so  that  the  mean  propeller  loading  remains  constant. 
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The  wakefields  vary  only  in  their  degree  of  circumferential  non- 
uniformity . 

The  propeller  blade  cavitates  mainly  in  the  vertical  upright 
position  because  the  low  inflow  velocities  present  there  cause 
the  blade  loading  to  increase.  Hence  the  inflow  velocity  in  this 
region  is  a major  factor  influencing  the  amount  of  cavitation  on 
the  blade,  and  an  appropirate  measure  of  the  wakefield  non-uniform- 
ity (as  it  influences  cavitation)  is 

Minimum  Velocity  at  0.9  Radius  (At  6=0°) 
Circumferential  Average  Velocity  At  0.9  Radius  5 


denoted 


^VMIN/VAVG^ '9R  ' 


If  this  quantity  is  unity  then  the  wake  will  be  perfectly  uniform 
and  the  transient  cavitation  will  go  away  completely.  As  this 
quantity  becomes  smaller,  the  intermittent  blade  loading  will  in- 
crease, which  will  cause  the  volume  of  unsteady  cavitation  to  in- 
crease . 

Figure  6.17  shows  the  estimated  cavitation  source  levels  for 
the  baseline  propeller  running  in  different  wakefields.  The  cor- 
relation between  ^mj^^AVG^  an^  sour>ce  level  is  clear.  This 
confirms  our  previous  conclusion  that  wakefield  non-uniformities 
play  a primary  role  in  determining  the  cavitation  source  strength 
of  merchant  ship  propellers. 
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7.0  COMPARISON  OF  PREDICTED  AND  MEASURED  ACOUSTIC  SIGNATURES 

Comparisons  between  measured  and  predicted  cavitation  source 
strengths  are  given  here  for  two  large  tankers.  The  acoustic 
measurements  were  made  by  Cybulski  (Ref.  18).  The  data  for  these 
ships  is  rather  sketchy,  and  we  have  had  to  estimate  the  wakefield 
and  some  of  the  propeller  design  specifications. 

Since  the  WORLD  DIGNITY  and  the  LAGENA  are  both  approximately 
the  same  size  and  shape  we  have  assumed  the  wakefields  to  be  the 
same,  for  lack  of  any  more  detailed  information.  The  wakefield 
was  estimated  from  Ref.  17,  and  is  the  same  as  the  wakefield  used 
for  the  sensitivity  analysis  (Pigs.  6.1  and  6.2). 

7.1  WORLD  DIGNITY 

The  available  data  for  WORLD  DIGNITY  is  listed  below. 

Length  1100  ft 

Beam  175  ft 

Draft  68  ft 

Tonnage  271,000  tons  deadweight 

Horsepower  ~30,000 

Maximum  Speed  16.0  knots 

Prop  Diameter  30  ft 

Prop  Depth  to  g ~50  ft 

No.  of  Blades  5 

Maximum  RPM  86 

Our  estimates  of  some  of  the  other  propeller  design  parameters 
are 

Pitch  Distribution  - reduced  pitch  at  blade  tip  relative  to 

optimum  pitch  distribution 

Expanded  Area  Ratio  - .60 

Percent  Lift  by  Camber  - 95$ 


w 
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The  predicted  source  level  for  the  first  three  harmonics  of  blade 
rate  are  shown  in  Fig.  7.1,  along  with  the  source  levels  as  infer- 
red from  acoustic  measurements  (Ref.  18).  We  have  not  shown  pre- 
dicted levels  at  the  higher  harmonics  because  the  influence  of 
bubble  dynamics  and  the  tip  vortex  cavity  volume  on  the  higher 
harmonics  have  not  yet  been  investigated.  We  feel  that  a compar- 
ison between  predicted  and  measured  levels  at  these  higher  harmon- 
ics would  be  premature  at  the  present  time. 

It  is  seen  that  the  predicted  blade  rate  level  is  12  dB  below 
the  measured  level,  while  the  second  and  third  harmonics  are  in 
fairly  good  agreement  (1-2  dB  too  high).  The  discrepancy  at  blade 
rate  may  be  due  to  hull  vibration  augmenting  the  radiation,  as 
discussed  in  Section  7.3- 

7.2  LAGENA 

The  major  characteristics  of  the  tanker  LAGENA  are  very  sim- 
ilar to  those  of  the  WORLD  DIGNITY,  and  hence  will  not  be  repeated 
here.  The  major  difference  is  that  the  LAGNEA  has  a 4-bladed  prop, 
while  the  WORLD  DIGNITY  has  a 5-bladed  propeller.  This  will  af- 
fect the  blade  rate  frequency. 

The  comparison  between  predicted  and  measured  source  levels 
for  the  LAGENA  is  shown  in  Figure  7.2.  Again  the  discrepancy  is 
largest  at  blade  rate  (17  dB)  and  decreases  with  increasing  har- 
monic number. 

7.3  Influence  of  Hull  Vibration 

There  are  several  factors  which  suggest  that  hull  vibration 
of  very  large  ships  is  augmenting  the  radiation  from  the  cavity 
volume  fluctuations  at  low  frequencies  (blade  rate  and  twice  blade 
rate ) : 
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o The  amount  of  cavitation  necessary  to  cause  the  measured 
blade  rate  level  is  shown  in  Figs.  7.1  and  7.2.  For  both 
ships,  this  amount  of  cavitation  is  far  above  that  usually 
observed  on  full  scale  ships.  Standard  naval  architect- 
ural practice  would  also  indicate  that  this  amount  of 
cavitation  is  excessive  and  would  probably  lead  to  blade 
erosion  and  possibly  thrust  breakdown. 

o Since  the  blade  rate  frequency  does  not  go  down  with  in- 
creasing ship  length  fast  enough  to  keep  fL  = constant, 
where  f = blade  rate  frequency  and  L = ship  length,  the 
blade  rate  frequency  (on  a non-dimensional  basis)  is  high- 
er for  very  large  ships  than  for  smaller  ships.  As  the 
frequency  increases  the  modal  density  for  most  structures 
will  increase.  Thus  the  probability  that  an  efficiently 
radiating  vibrational  mode  will  occur  at  blade  rate  fre- 
quency is  larger  for  big  ships. 

It  is  suggested  that  the  influence  of  hull  vibration  on  rad- 
iation be  investigated  before  exercising  the  cavitation  prediction 
model  to  predict  ocean  ambient  noise  due  to  propeller  cavitation. 
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8.0  CONCLUSIONS 

a)  The  present  method  of  calculating  sheet  cavitation  volume- 
time histories  appears  to  be  adequate  for  determining  the 
cavitation  source  level  at  the  first  few  harmonics  of 
blade  rate. 

b)  The  major  cause  of  cavitation  noise  at  the  higher  harmon- 
ics of  blade  rate  is  the  highly  non-linear  relationship 
between  blade  loading  and  cavitation  volume. 

c)  The  details  of  the  volume  history  will  affect  the  genera- 
tion of  noise  at  the  higher  harmonics  of  blade  rate.  The 
influence  of  bubble  dynamics  and  the  volume  associated 
with  a cavitated  tip  vortex  will  have  to  be  investigated 
before  we  can  confidently  predict  the  source  level  at 
higher  harmonics  of  blade  rate. 

d)  The  major  design  variables  which  influence  the  cavitation 
source  level  are: 

o diameter 

o expanded  area  ratio  (EAR) 
o percentage  lift  by  camber 

The  major  operational  variables  which  influence  cavita- 
tion source  level  are: 

o wakefield  non-uniformities 

o ship  speed  (as  a percentage  of  design  speed) 
o shaft  immersion 

e)  The  estimation  of  the  wake  non-uniformities  which  the 
full-scale  propeller  experiences  is  extremely  difficult. 
Although  wake  surveys  are  commonly  made  on  ship  models, 
the  full  scale  effective  wakefield  may  be  quite  different 
from  that  measured  during  a model  wake  survey  for  two 


reasons : 
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o The  large  Reynolds  number  difference  between  model  and 
full  scale 

o The  pressure  field  associated  with  a working  propeller 
substantially  alters  the  flow  over  the  stern  of  the  ship, 
and  hence  the  non-uniformities  of  the  wakefield. 

f)  The  presence  of  the  hull  may  affect  the  radiation  from 
the  cavitation  volume  fluctuations.  The  effect  may  be 
either  to  augment  or  decrease  the  radiation,  depending 
on  the  frequency  and  the  construction  of  the  ship.  This 
may  explain  the  extremely  high  blade  rate  source  levels 
observed  on  large  tankers. 

g)  All  of  the  propellers  investigated  in  the  present  work 
showed  significant  cavitation  only  when  the  blade  was  near 
the  vertical  upright  (0=0°)  position.  This  implies  that 
the  proper  source  depth  for  cavitation  volume  fluctuations 
is  near  thetip  of  the  blade  in  the  upright  position, 
rather  than  at  the  shaft  centerline. 
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9.0  RECOMMENDATIONS 

a)  The  Influence  of  bubble  dynamics  and  tip  vortex  cavita- 
tion on  the  cavitation  noise  at  higher  harmonics  of  blade 
rate  should  be  investigated. 

b)  The  effect  of  hull  vibration  on  the  radiation  from  a 
cavitation  volume  source  must  be  determined,  either 
analytically  or  experimentally. 

c)  The  numerical  discretization  errors  in  the  present  com- 
puter program  which  tend  to  generate  high  harmonics  need 
to  be  remedied. 

d)  Any  improvements  to  the  computer  program  which  will  re- 
duce the  cost  of  using  it  should  be  incorporated  before 
exercising  the  program  to  predict  the  ensemble  character- 
istics of  the  ocean  ambient  noise  due  to  propeller  cavi- 
tation. 
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